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Abstract 

We establish both upper and lower bounds on the electric field in the case 
where two circular conductivity inclusions are very close but not touching. We 
also obtain such bounds when a circular inclusion is very close to the bound- 
ary of a circular domain which contains the inclusion. The novelty of these 
estimates, which improve and make complete our earlier results in is that 
they give an optimal information about the blow-up of the electric field as the 
conductivities of the inclusions degenerate. 
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1 Introduction and statements of results 

The purpose of this paper is to set out optimal gradient estimates for solutions to the 
isotropic conductivity problem in the presence of adjacent conductivity inclusions 
as the distance between the inclusions goes to zero and their conductivities degen- 
erate. This difficult question arises in the study of composite media. Frequently in 
composites, the inclusions are very closely spaced and may even touch; see [S]- It 
is quite important from a practical point of view to know whether the electric field 
(the gradient of the potential) can be arbitrarily large as the inclusions get closer to 
each other or to the boundary of the background medium. 

There have been some important works on the estimates of the gradient of the 
solution to the conductivity problem in the presence of inclusions. For finite and 
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strictly positive conductivities, it was shown by Bonnetier and Vogelius in |1U| that 
the gradient of u remains bounded for circular touching inclusions of comparable 
radii. Li and Vogelius showed in ^1 that Vu is bounded independently of the 
distance between the inclusions B\ and B2, provided that the conductivities stay 
away from and +00. It is worth mentioning that the result of is much more 
general: it holds for arbitrary number of inclusions with arbitrary shape. This 
result has been recently extended to elliptic systems by Li and Nirenberg in [B] , 
On the other hand, for two identical perfectly conducting circular inclusions (with 
k\ = ki = +00) which are e apart, it has been shown in |B] (see also JZj and [T2"] ) 
that the gradient generally becomes unbounded as the distance e approaches zero. 
The rate at which this gradient becomes unbounded has actually been calculated in 

for a special solution. For this special solution, the rate turns out to be er 1 / 2 . 
In H3, a lower bound for the gradient of the solution to the conductivity problem 
for arbitrary conductivities, possibly degenerating, has been obtained. One of our 
objectives in this paper is to improve this bound and complete it by deriving an 
optimal upper one as well. 

In this paper we consider the following two situations: when two circular con- 
ductivity inclusions are very close but not touching and when a circular inclusion is 
very close to the boundary of the domain where the inclusion is contained. These 
two simple two-dimensional models illustrate very well the feature of our estimates. 
We believe that they extend to arbitrary-shaped inclusions if their contact reduces 
to a point. 

To describe the first situation we consider in this paper, let B\ and B2 be two 
circular inclusions contained in a matrix which we assume to be the free space M. 2 . 
For i = 1,2, we suppose that the conductivity ki of the inclusion Bi is a constant 
different from the constant conductivity of the matrix, which is assumed to be 1 
for convenience. The conductivity ki of the inclusion may be or +00. The zero 
conductivity indicates that the inclusion is an insulated inclusion while the infinite 
conductivity indicates a perfect conductor. We are especially interested in the case 
of extreme conductivities ki — ► +00 or ki — > 0. 

Given an entire harmonic function H, the first conductivity problem we consider 
in this paper is the following: 

Vfl + V (ki - l) X (Bi))Vu = in E 2 , 

< i^,2 ' (1) 

u(X) -H(X) = 0{\X\~ l ) as|X|^+oo. 

The electric field is given by Vu, where u is the solution to Q and represents the 
perturbation of the field VH in the presence of the two inclusions B\ and B2. For 
applications to the theory of composite materials, it is particularly important to 
consider the case when V-ff is a uniform field, i.e., H(X) = A ■ X for some constant 
vector A. The equation Q can be rewritten in the following form to emphasize the 
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transmission conditions on dBj, i = 1,2 



An = in n \ (dBi U dB 2 ) 

u\ + = u|_ on dBi, i = 1, 2, 



<9w 

Of 



on d-Bj, i = 1,2, 



du 
du 

u{X) - H{X) = 0{\X\~ l ) as|X|^+oo. 

Here and throughout this paper the subscript ± indicates the limit from outside and 
inside the domain, respectively. If ki = 0, then the transmission condition on the 
normal derivatives of u should be replaced with ^|+ = on dBi, while if ki = +oo, 
it should be replaced with u = constant on Bi. 

As has already been said, we are interested in the behavior of the gradient of the 
solution to the equation Q as the distance between the inclusions B\ and B 2 goes 
to zero for arbitrary conductivities k\ and k 2 , possibly degenerating. 

To state the first main result of this paper, let us first fix notation. For i = 1,2, 
let Bi = B(Zi,ri), the disk centered at Zi and of radius rj. Let Ri, i = 1,2, be the 
reflection with respect to dBi, i.e., 

Ri(x) ■■= r h x S Z tf + Zi > i = 1 ' 2. 

It is easy to see that the combined reflection R\R 2 and R 2 R\ have unique fixed 
points. Let / be the line segment between these two fixed points. Let Xj, j = 1, 2, 
be the point on dBj closest to the other disk. We also let 



min(ri,r 2 ), r max := max(ri,r 2 ), r* := a/ (2r 1 r 2 )/(ri + r 2 ), 

h + 1 , 1 

A,- := —r, 7T, i=l,2 and r : = 



2(ki-iy ' ' 4A!A 2 ' 

We obtain the following result on the blow-up of the gradient. 

Theorem 1.1 Let e := dist(B\,B 2 ) and let iM^ and j = 1,2, be the unit 

normal and tangential vector fields to dBj, respectively. Let u be the solution of 0). 

(i) If e is sufficiently small, there is a constant C\ independent of k\, k 2 , r\, r 2 , 
and e such that 

C 1 mf\(VH(X),u^(X j ))\ 

<|Vu| + (A^)|, j = l,2, (2) 



1 - t + (r*/r m i n )y/e 
provided that k\, k 2 > 1, and 

C 1 ml\(VH(X),T^(X,))\ 



Vt/ < |V«|+(^»U i = 1, 2, (3) 



1 - t + (r*/r min )^/e 
provided that k\, k 2 < 1. 



3 



(ii) Let Q be a bounded set containing B\ and Bi- Then there is a constant C2 
independent of k\, ki, T\, r<i, e, and Q such that 



C2\\VH\\ L oo(Q-) 

1 - |r| + (r*/r max )Ve' 



|V«|U«» ( n) < q 1,, , jz /„ ~ 7r - ( 4 ) 



Note that if H(X) = A ■ X for some constant vector A, which is the most 
interesting case, then the quantity 

(VH(X),^\X 1 )) = (A,^(X J )), 

and hence it does not vanish if we choose A appropriately. 

Theorem 11.11 quantifies the behavior of Vu in terms of the conductivities of the 
inclusions, their radii, and the distance between them. For example, if k\ and k 2 
degenerate to +00 or zero, then r = 1 and hence © and (|3J) read 

1 -£^<V«(X i )l J = 1,2, |V«|U- m < (5) 

for some positive constants C{ and C' 2 , which shows that Vu blows up at the rate 
of e -1 / 2 as the inclusions get closer. It further shows that the gradient blows up at 
X\ and X2, Xj for j = 1,2, being the point on dBj closest to the other disk. 

The lower bounds in (J2J) and (jSJ) are improved versions of the one obtained in [Hj. 
The proofs of our new estimates make use of quite explicit but nontrivial expansion 
formulae, originally derived in -5]. They are achieved by using a significantly different 
method from ^0], PS]- an d |Sj- 

Another interesting situation is when the inclusion is very close to the boundary. 
To describe this second situation, suppose that 0, which is a disk of radius p, contains 
an inclusion B, which is a disk of radius r. Suppose also that the conductivity of 
is 1 and that of B is k 7^ 1. The conductivity problems considered in this case are 
the following Dirichlet and Neumann problems: for a given / e C 1,Q (<9f2), a > 0, 



V • (1 + (fc - l)x(B))Vu = in Q, 
u = f on dO,, 



(6) 



and for a given g £ C a {dVt) 



V • (1 + (k - l)x(B))Vu = in O, 

— = g on dil . 

ov 

To ensure existence and uniqueness of a solution to 0, we suppose that j dQ g = 
and Jqq u = 0. Let X\ be the point on dB closest to d£l and X2 be the point on 
80, closest to dB, and let Rb and i?Q are reflections with respect to dB and dd, 
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respectively. Let Pi and P2 be fixed points of RbRq, and R^Rb, respectively, and 
let Ji be the line segment between P\ and X\ and J2 that between P2 and X2. 

The second main result of this paper is the following triplet of estimates for the 
gradient of the solutions to © and (J7J). Let X>n(/) and Sn(g) denote the double 
and single layer potentials whose definitions are given in Section 2. For the Dirichlet 
problem we have the following theorem. 

Theorem 1.2 Let 

k—1 



e := dist(B, dQ), a 



k + 1 V P r 

and let u be the solution to J^l . 

(i) If k > 1, then there exists a constant C\ independent of k, r, e, and f such 
that for e small enough, 

d inf \(VV n (f)(X),MXi))\ 

e 1 , T < |Vd+(Xi)|, (8) 

and 

Ci inf l(VC(,(/)(Jf),^(X 2 ))| 

!- J + fr .Vi <|V«I-(^)|. (9) 

Here vb and vq, denote the outward unit normal to dB and dVL. 

(ii) For any k ^ 1, there exists a constant C2 independent of k, r, and e such that 
for e small enough, 



livtill^m) < f 2|m , lcl '" ( T (10) 

For the Neumann problem the following theorem holds. 

Theorem 1.3 Let e, a, r* be defined as in Theorem M.IA 

(i) If k < 1, then there exists a constant C\ independent of k, r, e, and g such 
that for e small enough, 

d inf \(VSn(g)(X),T B (X 1 )}\ 

< iViil+TOI, (11) 



1 + cr + 4r*^e 
and 

Ci inf |(V5a(5)(X),T n (X 2 )}| 

e 2 ■= <|Vu|_(X 2 )|. (12) 

Here Tb and Th denote the positively oriented unit tangent vector fields on dB 
and dQ, respectively. 
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(ii) For any k ^ 1, there exists a constant C 2 independent of k, r, and e such that 
for e small enough, 



HV«|U~ ( n) < C2l , lg ! lc " (8 ^ - (13) 
" UL (U) ~ 1- \a\ +r*yfe y ' 

If Z is the center of £1 and if f(X) = A ■ X for some constant vector A, then 
T> n (f)(X) = \A ■ X for X e Q and V n (f)(X) = -^0^ for X G R 2 \ H, and 
hence we can achieve 

(VPn(/)P0, MXi)) + and (VV n (f)(X), v n (X 2 )) ± for any X, 

by choosing A appropriately. Likewise, if g := A ■ v on dft, then S n ( ff ) = -\A ■ X + 
constant. 

Theorem II .21 shows that in the case of the Dirichlet problem, if the inclusion is 
a perfect conductor (k = +00 and hence a = 1), then 

C\ „ a 

r*ye v ' r*ye 

for some positive constants C{ and C^. Thus Vu blows up at the rate of e _1//2 as 
long as the magnitude of r is much larger than that of e. It also shows that the 
gradient blows up at the points X\ and X 2 - On the other hand, for the Neumann 
problem, according to Theorem 11.31 the situation is reversed: Vu blows up for an 
insulator. If r is of the same order as e, then r* ~ -K= and hence Vu does not blow 
up. In fact, it stays bounded and an asymptotic expansion of the solution as e — ► 
can be derived. See for instance ^ HJ [7j for this. 

One may think that Theorems 11.21 and 11.31 can be easily derived from Theorem 
11.11 by reflection or conformal mapping. This is, as it will be shown, far from being 
true. The proofs of Theorems 11.21 and 11.31 require delicate analysis and careful and 
tricky estimates. 

In this paper we only deal with the two dimensional case. It seems challenging to 
obtain similar results in three dimensions. At this moment it is even not clear what 
the blow-up rate of the gradient would be in three dimensions. In this direction, we 
have recently established in 2: that, unlike the two-dimensional case, if the inclusions 
are grounded conductors (the Dirichlet boundary condition on the inclusions is set 
to be zero) and of spherical shape, the gradient stays, to our surprise, bounded 
regardless of the separation distance between them. 

This paper is organized as follows. In Section 2 we review some preliminary facts 
on layer potentials and representations of the solution to the conductivity problem 
obtained in jSj. Theorem ll.il is proved in Section 3, and Theorems 11.21 and 11.31 in 
Section 4. Although our results hold for special cases, we believe that they extend 
to arbitrary-shaped inclusions if their contact reduces to a point. 
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2 Preliminaries 

To make our paper self-contained and our exposition clear, we introduce our main 
tools for studying the conductivity problems and collect some preliminary results 
regarding layer potentials. We also linger over a description of some results from our 
earlier papers. The material in Lemmas l2.5l and l2.6l is however, up to our knowledge, 
new. 

Let D be a bounded Lipschitz domain in M? and Sjjcf) and £>d^> denote the single 
and double layer potentials of a function <p G L 2 (dD), namely, 



if the limits exist. Here and throughout this paper vx is the outward normal to dD 
at X and ( ,) denotes the scalar product in M 2 . The following jump relations for 
the single and double layer potentials are well-known 118| . 

Lemma 2.1 For (ft G L 2 (dD) we have 




For a function u defined on R 2 \ dD, we set 





(14) 



(v D 4>)\ ± (x) 



(T-+lC D )<f>(X) 



a.e. X G dD, 



(15) 



where K,r> is defined by 




and K,* n is the L 2 - 



adjoint of ICd, i.e., 




Here p. v. denotes the Cauchy principal value. 



Note that if D is a two-dimensional disk with radius r, then 



(X - Y, v x ) 
\X-Y\ 2 



— , VX, Y G dD, X^Y, 
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and hence 

JC D (f>(X) = K D (p{X) = -*- I <f>(Y)du, X G dD. (16) 

47rr ,/ SD 

When D is a disk, we define the function i?£>(/) for a function / by 

R D {f){X) :=f(R D (X)). 

Using the jump formula for the single layer potential, we get the following lemma 
0- " 

Lemma 2.2 Let D be a disk in M? and let Rp denote the reflection with respect to 
dD. If v is a harmonic in D and continuous on D then 



S D 



fdv_ 



1 



(R D v)(X) + C, XgR 2 \D, 



(17) 



where C is some constant. Analogously, if v is harmonic inR 2 \ D, continuous on 
R 2 \ D, and v(X) -> as \X\ -> 0, then 



S D 



0v 
dv 



dD 



(X) = ^(R D v)(X) + C, X £ D, 



(18) 



for some constant C. 

We need the following lemma which was first proved in jSJ. 

Lemma 2.3 Suppose that B\ and B2 are two disjoint disks and let Ri = R Bi be the 
reflection with respect to dBi. Then the solution to Q) is represented as 



U {x) = H(x)+s Bl Mx)+s B2 ^ 2 (x), xen, 



(19) 



where ipi G L^{dBi), i = 1,2, is the unique solution to the system of integral equa 
Hons 



Q v (l) IdBi Q v \8B. 



— — L R on dBi, 1 = 1,2, i^l, 



(20) 



with Xi = 2{fcT~n an d v ® ^ s the outward unit normal to dBi . Moreover, the potentials 
if 1 and if2 are explicitly given by 



tp 2 



1 + °° 1 

- V - 







Ai ^ n (4AiA 2 ) m duO) 

771=0 

i_ ^ 1 a 

A 2 ^ n (4AiA 2 ) m dvW 

771=0 



{R2R1TV ~ W R2)H 



{RiR2) m (I- w Ri)H 



8B 2 



(21) 



where the series in \21\) converge absolutely and uniformly. 
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In Lemma 12 .31 the space L^(dBi) denotes the set of all g G L 2 (dB{) having mean 
value zero: f dB g = 0. The following lemma from jS] is also of use to us. 

Lemma 2.4 Let u be the solution of 0j and let H be a harmonic conjugate to H . 
Let v be the solution to the conductivity problem 



Then 



' v{i+£(±-i)x(i*))v 

i=l,2 1 
v(X)-H(X) = 0(\X\- 1 ). 

dv 



in 



(22) 



du 
df 



on dB{ , i = 1, 2. 



We now turn our attention to the second situation, i.e. problems (jHJ) and 0, 
when both Q and B are disks. We first note that T>q/ is C ,a in f2 and M? \ since 
/ £ C 1,Q (<9il). It was shown in [LSI I14j that the solution u to the problem Q for a 
fixed Dirichlet data / is given by 

du , 



(23) 



where (/? with mean value zero satisfies the integral equation 

d 



on 



OB, 



with A 



(\I-K.* B )<p = j-(V n (f)-Sn(g)) 
k + \y Since B is a disk, it follows from (HI that /C^ = on L%(dB) 



2{k - 1) ' 



and hence 



On the other hand, g 



d_ 

du 



{Vn(f) - S n {g)) ondB. 



(24) 



f^lan yields 
(9 



Since ^S n (g)\- 



g = ^(V Q (f)-S n {g)+S B (<p))\_ on dU. 
■ {-\L + KL* n )g and $7 is a disk, J^Sn(g)\- = -\g on 95. Thus 



we get 



\g = |;(%(/) + Sb(vO)|_ on ^. 



(25) 



It then follows from (|24j) and (|25() that 5 and 9? are the solution of the following 
system of integral equations 

1 d{s B ip) d(v n f) 
2 9 



d(Sng) 
du B 



du n 

d(Vnf) 
dv B 



on d£l, 
on dB. 



(26) 
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Observe the similarity of (|26|) to H2U[) , Using the same argument as the one intro- 
duced in deriving (|21|) . one can show that the following lemma holds. 

Lemma 2.5 Let g and <p be the functions given in h2°J\) . Then g and (p are given 
by 

a = 2 T,j^^ r b^T(i- Tx RbW\ on on, 

1 1 ft 

v = x^j2xrd^ [iRnRBr{I ~ Rn)I)nf] ondB ' 

in=0 

where the series in \21\) converge absolutely and uniformly. 

Proof of Lemma 12. 51 The convergence of the formula (|27|) will be proved in the 
course of proving Theorem 11.21 in Section 4. 

We first prove that for (hi, h 2 ) G L^(d£l) x Lq(8B) there exists a unique solution 
(g,cp) G Ll(dQ) x Ll(dB) such that 

1 d(S B <p) 

2 9 dv^~ = ° n 

> . d(S n g) 

\ip H = I12 on oB. 

dv B 

Since B is away from dQ., the operator (g,<p) — * (— , 9 ^g^ ) is a compact 

operator on Lq(cK2) x Lq(8B). Thus, by the Predholm alternative, it suffices to 
prove that if (hi,h 2 ) = (0,0), then the solution (g,<p) = (0,0). In order to do this, 
suppose that 

1 d(s B y) 

-g = on dQ, 

2 dm 

Xip H = on oB. 

Then the function u = Sn(g)+S B (ip) in f2 is a solution of V-(l+(jr — i)x(B))Vu = 
in f2 and satisfies = on dil. This implies that u is constant in f2, and hence 
S B (cp) is harmonic in f2. It then follows from the jump formula (|14[) that 99 = and 
therefore, Sn(g) = constant in $7. Hence, 5 = — 2 8 g" g |_ = 0. 

We now prove that the pair (5, <£>) given by ((2"T|) satisfies (|26j). Observe that the 
function 

(R B Rn) m (L-^R B )(V Q f)(X) = (V Q f)((R n R B ) m ^ 
is harmonic in M? \ B and approaches to 

(V n f)((R n R B ) m - l R n (Z)) - ^(V n f)(R B (R n R B ) m - l Rn(Z)) 
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as \X\ — > +00, where Z is the center of B. Since 5fj(l) is constant in Q, it follows 
from (JT7J) that 

m=0 
m=0 

Likewise one can show that 

m=0 

Thus (5, (/?) satisfies (|2f))) and the proof is complete. □ 

The representation of g and 93 given in (|27|) can be simplified using the relation 
RnV n f(X) = V n f(R n (X)) = -V n f(X) + constant, l£K 2 \ (28) 
which follows from (|15|) and Q16[) since SI is a disk. Using l|28[). we then compute 

—\(F! D F!^ m (T - 

2A 



m=0 



m=0 



1 


d 


(2A) m 


dm 


1 


d 


(2A) m 




1 


d 


(2A) m 





<Ep^il(«A)"a/] + ^%/ on afi. (29) 

Likewise we can show that 

2 + °° 1 «9 

m=0 

The following lemma is also of importance to us. 

Lemma 2.6 Let u be the solution of Q) for g £ C a (d£l) and let G be the function 
satisfying §^ = g on d£l and J dn G = 0. Define v to be the solution of the following 
conductivity problem 



V-(l + (--l)x(-B))vu = inn, 



(31) 



v = G on d£l. 
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Then 



du 
df 



dv 
du 



on dB. 

Moreover, T>q(v\q^) is a harmonic conjugate to Sag in f2. 



(32) 



Proof. Let w be a harmonic conjugate of u in 0, \ B and B. Such a conjugate 
function exists in $7 \ B since J c ^da = for any simple closed curve C in f2 \ B. 
Moreover, since u is C 1,a , so is w. 
Define v by 



«P0 := 



(w(X), 
j &w(X) 



1951 



w da, X £ 5. 



(33) 



Then one can see from the Cauchy-Riemann equation and the transmission condi- 
tions on u that 



dv 
df 



du 
du 



= k 



du 
du 



dv 
df 



dv 


du 


du 


1 dv 


du 


+ = ~df 


+ = ~df 


k du 



Thus v defined by (|33j) is the unique solution to l[22j). and hence (|32|) holds. 
It follows from (fT3)l that 



and hence 



9(2?n(w|«i)) 



-v on 



an, 



1 du _ 1 

2 9T ~~ 2 9 



d(S n g) 



on 9f2. 



2dT 2 y 9^ 

Therefore Pn^lan) is a harmonic conjugate of Sng in IX This completes the proof. 



□ 



3 Proof of Theorem 11.11 

At this point we have all the necessary ingredients to prove Theorem 11.11 As has 
been said, the lower bound in (J2J) and (JHJ) is an improvement of the one obtained in 

Recall that there are two disks Bj = B(Zj,rj), j = 1, 2, inside f2, and that Rj is 
the reflection with respect to dBj. We suppose that both centers Z\ and Z^ are on 
the x-axis. 
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If X is on the x— axis, i.e., X = (x,0), straightforward calculations show that 

DR i (X)=g l (X) 



-1 
1 



and 



where 



Therefore, 



V(R i f)(X) = Vf(R i (X))g l (X) 
9i{X) 



-1 
1 



\x-ztf 



, * = 1,2. 



V((R 2 Ri) m H)(X) 



2m 



l[g li (R li _ 1 ---R h (X)) 



i=l 



VH((i? lJ R 2 ) m (X)) 



1 
1 



(34) 
(35) 

(36) 



and 



V{{R 2 R 1 ) m R 2 H){X) 



92 {{RiR 2 ) m {X)) 



2m 



Hg k (R k _ 1 ---R h (X)) 



i=i 



VH(R 2 (R 1 R 2 ) m (X)) 



-1 N 
1) ' 



(37) 



times 



where (h, . . . , l 2m ) = (2, 1, 2, 1, . . . , 2, 1). 

Let u be the solution to the equation (JTJ). Combining (JT1J), (UHJ), (HHJ), and (|2"U1) 
yields 



<9w 



+ 



± 



(Aj ± -)v?j on &Bj, i = 1,2. 



Consequently, 



|Vu|±(AQ)| > 



du 



> 



Ai ± - 
2 



(38) 



Suppose that fci > 1 and A; 2 > 1. By (|2T|) . (|36|) and (|37jl. we obtain the following 
inequality: 

" Ai (4AxA 2 )- V 2A 2 J i 1 1 



> 



m=0 
1 1 -| 2_ 



Ai l - . a \ i 



^-inf |Vi? -i/Wl 



a 2 
4AiA 2 
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where a := (1 + 2(r*/r m i n )-y/e) . We then get from (|35|) that 

iv " Mxi)ia i- T+( , ( :/ rm ,„)^ y |vg - y '" 1 - 

Likewise we obtain 

|Vn| + (X 2 )| > —mi\VH-u^\. 

Thus the lower bound in © is now derived. 

If both k\ and fc 2 are less than 1, then one can use Lemma |2.4I to obtain the 
lower bound in (J3J). See jH] for the details. 

We now prove (JIJ. We need the following lemmas. 

Lemma 3.1 Let r max , r m i n , r*, and e be as in Theorem M.ll If e is small enough, 
then for any X £ B\ and n > 8r^/y/e we have 

\(RiR 2 ) n {X) - Z 2 \ > r 2 (l + -^V~e), 

and 

\{R 2 R x ) n R 2 {X) - Z X \ > n{l + ^y/e). 

"max 

For any X E B 2 and n > Sr^/y^ we /icwe 

|(i? 2 i? 1 ) n (X)-Z 1 |>r 1 (l + -^^), 

and 

KR^TR^X) - Z 2 \ > r 2 (l + -^-y^). 

"max 

Proof. After a translation and a rotation if necessary, we may assume that B x = 
B((0,0),n) and B 2 = B{(n + r 2 + e,0),r 2 ), i.e., Z x = (0,0) and Z 2 = (n + r 2 + e,0). 
It is easy to show that the fixed points of the combined reflections R 2 R X and R X R 2 
are the points (xj,0), for i = 1,2, where x x and x 2 are the roots of the quadratic 
equation 

n+r 2 + e^j x 2 + —r\ — (n + r 2 + e) 2 ^ x + r\ {r x + r 2 + =0. 
Then, as e goes to zero, 



xi = r 1 - A /^ 2 - v ^ + 0(6), x 2 = n+ J*p-yTe + 0(e), 
w n + r 2 V H + ^2 
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and 



-r*Ve < |n - xj\ < 2r*^/e, j = 1,2. 
Let Xi = (n,0) and (t n ,0) = ( J Rii? 2 ) n (^i). We have 



(39) 



(n + r 2 + e) 



n + r 2 + e - t r , 



and hence 



|Wi - x i\ 



i 



(n + r 2 + e 
|*n - zi 



ri + r 2 + e - t n 

nr 2 



(n + r 2 + e) 



7 2 



n + r 2 + e - n 



(n + r 2 + e)(ri + r 2 + e - i n ) - 
rir 2 



< 



Thus, 



(n + r 2 + e)(n + r 2 + e — x\) - r\ 

\tn ~ Xi\ 
l + # ' 



|*n - xi\ < ifn>-^=. 



(40) 



Observe that | ( J Rii? 2 ) n (Xi) - Z 2 | = n - i n + r 2 + e. Therefore, it follows from JSU) 
and (gDJ) that 



\(R 1 R 2 ) n (X 1 )-Z 2 \>r 2 + ^ ifn>^. 

z v e 



(41) 



Let X E B\. For any < t < ri such that |X — Z 2 | > |(t, 0) — Z 2 \, we can easily 
see that 

\R 2 (X)-Zx\ > \R 2 (t,Q) -Z x \ and \RiR 2 (X) - Z 2 \ > |i?ii? 2 (t, 0) - Z 2 \. 

Since RiR 2 (t,0) = (s,0) for some s satisfying < s < r\, then, by repeating the 
above inequalities, we have for each positive integer m, 

\{R 1 R 2 ) m (X)-Z 2 \ > KRrf^^O) - Z 2 \, 
\R 2 (RiR 2 ) m (X) - Z X \ > l^^^r^O) -Z X \. 
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In particular, for the case t = n, we obtain that 

KR^riX) - Z 2 \ > \(R 1 R 2 ) m (X 1 ) - Z 2 \, 

\R 2 (RiR 2 ) m (X) - Z x \ > \R 2 {R 1 R 2 ) m {X 1 ) - Zi|, Vie B x . 

Similarly, we have 

\{R 2 R x ) m (X) - Z x \ > \(R 2 Ri) m (X 2 ) - Z X \, 

\Ri(R 2 Ri) m (X) - Z 2 \ > \R 1 (ReR 1 ) m (X 2 ) - Z 2 \, VIeB, 

Therefore the first two inequalities in Lemma lM . 1 1 follow from (|41jl. The second pair 
of inequalities can be derived by interchanging B\ and B 2 . This completes the proof. 
□ 



Lemma 3.2 (i) Let X\ = (ri,0). We have 

g2((RiR2) m (X 1 )) ) gxHR^TR^)) > ) V m e N. 

1 + 8(r*/r min )Ve 



(42) 



(ii) For all X € B\, we have 



g 2 {(R 1 R 2 ) m (X))<l, VmGN, 

1 (43) 

92 {(RiR 2 ) m (X)) < —— — , V m > 8n/V~e, 

i + ir*/r max j \/ e 

and similarly, for all X € B 2 , we have 

'g 1 ((R 2 R 1 ) m (X)) < 1, VmeN, 

1 (44) 

gi((R 2 Ri) m (X)) < — — — , V m > Sn/V~e. 

l + U*/ r maxJv e 

Proof. Since (RiR 2 ) m (Xi) is between X\ and X 2 where X\ = (xi,0) and X 2 = 
(:r 2 ,0) (x\ < x 2 ) are fixed points of R 2 R\ and RiR 2 , respectively, we have 



g^RxR.riX,)) = — -— 7 , > 



Kfli.Ra)" 1 ^!)-^! 2 " |X!-Z 2 |2 



> =■ > 



(n + r 2 + e - xi) 2 (r 2 + 2r*^/e + e) 2 1 + 8(r*/r m i n )y/e 
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The second inequality in (|42(1 can be proved in exactly the same way. Lemma 13,11 
and the definition of gi give the upper bounds 1)43)1 and (|44|) . □ 



To establish our upper bound we first observe that since u{X) — H(X) — » as 
|_X"| — > +00, |V(u — H)\ attains its maximum on either dB\ or 8B2, and hence 



< \\V(u - -H')|+||L°°(aBiu9B 2 ) + IIV-^IIl 00 !^) 

< \\^ u \ + \\l^ (dB 1 UdB 2 ) + 1 1 V-H" llioo^), 



(45) 



and 



1 1 Vlt Hloo^u^) < || Vli|— \\L°°(dB 1 UdB 2 )- 

We also have from (fT9j) and (|38|) that 



|Vm|±|| L oo( 9b .) < (| A*| + -) 11^11x00(9^) + 



du 



dT 



(46) 



Let N be the first integer larger than Sr^/y/e. It then follows from (|21j) . 1)4311 . 
and (@3I) that 

MX)| < ||ViT|Uoo( BlUB2) ^f J] — - ,_ (!+ ' 



^<7V 



|4At A^| m v * ' 2|A 2 



) 



+ |4AiA 2 r ^ |4A X A 2 ' ! ' 

m=0 



2m 



for any X G dB\, where b := 1/(1 + (r*/r max )y/e). Thus, for each X £ dB\, 



M*)l < 



C||VH||l°o(b 1 ub 2 ) 



1 



+ 



1 - It 



|8r„/v^ 



I Al| V 1 ~ l r l + r */ r max\/e l-| r l 

|V-ff||ioo( BlUB2 ) 



I Ax I (1 - |t| + r*/r max .y/e) ' 
for any for any X E <9-Bi . Similarly, we have the desired estimate for tp2 (X) , X £ <9£>2 



and hence, 



3u 



C||V-ff||2/x>( Bl uB 2 ) 



du L°°(dB 1 VJdB 2 ) ^ ~ l T l + r */ r max\/e 

To estimate du/dT we use Lemma 12.41 Let be a harmonic conjugate of H 
and v be the solution to 1)22)1 . Then by 1)47)1 . we have 



< 



L° a (dB 1 UdB 2 ) 



C , ||^ 7 ^IU° (BiUS 2 ) 

1 - |r| + r t /r mEK \/e' 
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Since ||Vi7||£oo( BlUB2 )||V-ff||£oc( BlUB2 ), it follows from Lemma that 
du 



dT 



< C \\ VH \\^^B 2) (48) 
L°°(dB 1 UdB 2 ) * l r l + r */ r maxVe 



Combining (fiT|) . and yields the upper bound in (@J) and completes 

the proof. 

4 Proofs of Theorems 11.21 and 11.31 

We suppose that Q = B((0,0),p) and B = B((p — r — e,0),r) after rotation and 
translation if necessary, and that e <C p — r. The conductivities of f2 and £? are 1 
and k, < k ^ 1 < +00, respectively. Let 

2 2 
g n (X) = j^, g B (X) 



|X| 2 ' 7 \X- (p-r-e,0)| 2 ' 

The functions and g B play the roles of g\ and 52 in the previous section. Let Rq 
and R B be the reflections with respect to and B. 

Lemma 4.1 Lei Xi = {p — e, 0), i/ie point on B closest to dVL. For any positive 
integer n and X G B, we have 

gB(Rn(RBRn) n (X))gn((RBRn) n (X)) <g B (Rn(RBRn) n (X 1 ))gn((RBRn) n (X 1 )). 

(49) 

Proof. For any X = (x,y), one can easily see that 

2 2 

g B (Ra(X)) gu (X) - 



(5fi&-(P -'-«)) +(#) X +y 
p 2 r 2 1 



Since p - r - e > P — e ) it immediately follows that 

ffB (i?n(X))^(X) < <? B (RniX^gniXt), VX e 5. 
If X satisfies \X — (p — r — e)| < \(t, 0) — (p — r — e)|, with t > p — r — e, then 

|i? B ife(X)-(p-r-e)| < | J R Bj Rn(t,0)-(p-r-e)|. 
Using this fact repeatedly, we have 

|(i? B i? n ) n (X) - (p - r - e)| < \(R B Rn) n (Xi) - (p - r - e)|. 
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By combining this inequality with (|5U|) . we obtain that, for X G B, 
9B(Rn(RBRn) n (X))gn((RBRn) n (X)) <g B (Rn(RBRn) n (X 1 ))gn((RBRn) n (X 1 )), 
which completes the proof. □ 

Recall that Pi and P2 are the fixed points of the combined reflections RbRq. 
and RqRb, respectively. Observe that P\ G B. If Pi = (xj,0) for i = 1,2, then 
Xi {x\ < X2) are the roots of the quadratic equation 

(p — r — e)x 2 + (r 2 — p 2 — (p — r — e) 2 )x + p 2 (p — r — e) =0. 

It then follows that 



*i=P-\ —V~e + 0{e), x 2 = p+ Jj^L^ + 0(e). (51) 
V p — r y p — r 

Moreover, 

^< \xj-p\ <2^,j = l,2, 

for e small enough. 

As a direct consequence of (j49j) and (|51|). we have the following lemma which 
plays a crucial role in deriving the lower bound (JHJ). 

Lemma 4.2 For each positive integer n, the following inequality holds: 

gB(Rn(RBRn) n (X 1 ))gn((RBRnnx 1 )) > — (52) 

Proof. Since Pi G B, we have 

5B(^(^^)"(X 1 ))^((P B i? n ) n (X 1 )) > <7B(^(PB^)"( J Pi))^((^^) n (i :, i)) 

> 9B{Rn(Px))gn(Pi). 

On the other hand, it follows from (|50|) and (j51[) that 

9 B (Rn(Pi))gn(Pi) = — — > 



(^-xi^-r-e)) 2 l + 4r*^' 
The proof is complete. □ 

The following lemma is also of use to us. 
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Lemma 4.3 For each positive integer n, similarly to \5ty) the following holds: 

1 



gn(RB(RnRB) n (X 2 ))gB((RnRB) n (X 2 )) > 



1 + 4r*^/e' 



(53) 



Proof. Since the proof of Lemma 14.31 is parallel to that of Lemma 14.21 we very 
briefly sketch it. As before, since P 2 £ K 2 \ Q, we can show that 

9B((R^RB) n (X2))gn{(RBRn) n RB(X 2 ))>gB({RnRB) n (P2))gn{(RBRn) n RB(P2^ 

= 9B(P2)gn(R B (P2)). 

Since, as one can see easily, Rq(Pi) = P 2 we have 

9B(P2)gn(R B (P2)) = 9B(Rn(Pi))gn(Pi), 
and hence (|53|) follows. □ 



Next, we need the following lemma to derive the upper bound. 
Lemma 4.4 For each positive integer n, 

gB(Rn(RBRn) n (X))gn((RBRn) n (X)) < 1, IeB, 
and for n > . } r , 

— 4r*ye' 



gB{Rn{R B Rn) n {X))gn{{RBRn) n {X)) < 



1 + r*^' 



X G B. 



(54) 



(55) 



Proof. The inequality ((5*1)) is obvious. We only have to prove (f""5~|) . Let (i ra ,0) := 
{R B Rn) n {Xi). Then we have 



b n+l — p2 



k~0>-r-e) 
Recall that Pi = (xi,0). We then have 

\tn+l ~ x i\ 



+ p — r - e. 



< 



^-{p-r-e) ^-(p 

\ L n -^l 



r — e 



p 2 — (p — r — e)xi 



pr 



p 2 -{p-r - e)t n 



p—r 



1 + 
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If n > 5^-, we get 



and therefore 



1 / pr _ 

2 y p — r 



''' <V " ' ' ^2{f^r)^ : - (56) 



By Lemma IP and (f5U|). 

2 2 

9B(Rci(RBRa) n {X))g n ((R B Rn) n (X)) < 



1 

< 



1 + r*^e' 

which completes the proof. □ 

Since RbRq(X) is in £> for any I £ !!, we immediately obtain the following 
corollary. 

Corollary 4.5 For eac/t positive integer n, 

gB(Rn(RBRn) n (X))gn((RBRn) n (X)) < 1, le!), 

and for n > 

gB(Rn(RBRn) n (X))gn((R B Ra) n (X)) < — 1-^ , leH. 

We are now ready to prove Theorem 11.21 and Theorem 11.31 
Proof of Theorem \1.2l Straightforward computations yield that for X\ = (p — e, 0) 

V((i2 n i2 B ) m 2W)(*i) 

m— 1 

= [J gBiRniRBRnnX^gnttRBRnnX^VVnfttRBRnnx^). (57) 

n=0 

Since Ra(Pi) = P2 and Rb{P2) = A, (-R_B-Ro) n (^i) lies in J\, the line segment 
between Pi and X\, for each n. We may assume 

inf \(VV Q f(X),u B (X 1 )}\^0, 

since otherwise the estimate (JHJ) is trivial. If e is small enough, then the length of 
J\ is small and hence we may further suppose that 

(VVnf(X)^B(Xi)) = ^pV) 
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has the same sign for all X € J\. It then follows from (|52j) and l)57|) that 



dv B 

m—l 



n=0 

>a m inf 



<9x 



(58) 



where a := (1 + 4r* A /e) 1 . 

Suppose that k > 1 and hence A > 0. It follows from that 



\ 2^ fOU 



A ^ (2A) m &/ B 

m=0 



[(RnRBrVnfKX, 



> 



m=0 



f a \ m 
A ^„ V2A 



inf 



dx 



(X) 



> 



Ca inf^eJi 



1 -a + Ar*Je 



since cr = Here C is a constant which is independent of fe, r, and e. 
By (|26j). we have 



du 



dv B 



and hence we obtain 



_ dV n f 




dSng 


± 9z/ B 


± dv B 


± 


dv B 



<9u 



(*1 



Cinf 



> 



± = (\±-)<p on 55, 



9(T>af) f X \ 



(59) 



1 - 0- + 4r*Ve 



(60) 



In order to prove ©, we now estimate from below g{X2) where X2 = (p, 0), the 
point on <9il which is the closest to dB. Elementary computations show that 



d 

dm 



[{R B Rn) m Vnf]{X 2 ) 



m—l 



[1 gn(R B (RuR B ) n (X 2 ))g B ((R n R B ) n (X 2 ))(VV n f((R n R B ) m (X 2 )), u n (X 2 )). 



n=0 
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Since (i2o-Rs) n (X2) lies in J 2 , the line between X 2 and P 2 , for each n, we have as 
before 



9(X 2 ) 



> 



j +00 
9X 



I2A ^ V2A 

m=0 



a \ m 



+ 2 



• inf 

xeJ 2 



dx 



(X) 



> 



1 - <7 + 4r* v / e 



Therefore, we get 



du 



dm 



(X 2 



Cinf 



> 



xeJ 2 



g(gn/) 



9x 



(X) 



I- a + 4r*Jl 



We now prove 1)10(1 . Let iV be the first integer such that iV > 
get from Lemma 14.41 that 

\V((R Q R B ) m V Q f)(X)\ 



4r*Ve' 



(61) 
We then 



< II 9B(Rn(R B Rn) n (X))g n ((R B R n ) n (X))\VV n f((R B R n ) n (X))\ 

n=0 

[|V2?n/||j 



< < 



for all m, 



um ~ N if m > N, 



lL°°(f2) 

J|V2>n/IL«<ri)&' 
where b : = (1 + r*i/e) _1 - It then follows that for all X £ B 

1 + °° 1 



1 1 +°o I 

/ v — v / 1 \ m I x -« / b \ m 
<2||VPo/IU ( o)(E(^) V^tl 

m<Af 11 11 m=0 1 1 

<C||VPn/[U«(n)( ^ + 



l - M 



1 - Icrl +r*v/e 



< 



C||VPn/|| L oo (n) 
1 — \a\ + r*J~e 



(62) 



By Lemma I2~2l and (J3IJ), |V5s<p(X)| = g i j(X)|V5 B (^(i? i3 (X))| for l£Sl\B, and 
hence (|62() holds for all X £ Q, i.e., 



11 . g||VP n /|| LOO(n) 

l|V5 flV ||^ (n) < j^p^ • 



(63) 
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Since u is harmonic in B, it follows from (|59|) that 
fdu \ 

u(X) = -2Sq — )(X) + constant(l - 2\)S B ip(X) + constant, X e B. 



dv 

Since 2A - 1 = 2/(fc - 1), © gives 



IV7 II ^ g||VP n /|| LO c (n) 

K ' \k — 1|(1 — |<r| + r*^) 



(64) 



By dHJ and (JTSJ), 



+oo ^ 

S n g = -2 ^ 72X^( R n R B) m Vnf - V a f + constant. 

m=l ^ ' 

By Corollary 14.51 and computations similar to those in (|62j). we obtain 

||V«S^|| LOO(II) < 1-k|+r ^ • (65) 
By combining (|23|) . (|62|). and (|65j). we arrive at 

nv7 ii / „l|VI>n/IU-(n) , v 

l|V«|U<» ( n) < C Y^]^[+7%7i' (66) 

Since maps C 1,Q (50) into itself, it follows from the maximum principle that 

||V2W|| ioo(n) < ||VX> n /||L»(«i) < C\\f\\ c i, a{dQ) , 

and hence we get (|10j) . The proof of Theorem 11.21 is now complete. □ 

Proof of Theorem \1.3[ To prove Theorem II .31 we use Lemma 12.61 Let G be the C 1,a 
function such that ^ = g on <9S7 and Jq^G = 0, and v be the solution to Q31|l . 
Since Pq(G) is a harmonic conjugate to 5^5 in Q and 1R 2 \ $7 by Lemma f2.61 we 
have 

{VSn{g){X),T B {X 1 )) = -(VP n (G)(X),^(X!)), 

and 

(V«S n (g)(X),T n (X 2 )) = -(VV Q (G)(X),va(X 2 )), IeK 2 \ Oft. 

Thus (dU and JT2J) follow from ©, ©, and (j32J). 

As one can see in the proof of Lemma 12.61 v is a harmonic conjugate to u in 

Q \ B and — v is a harmonic conjugate to u in 1?. Therefore by (1661) we get 

nv7 ii ^ iiv7 ii ^ C|l v '5n5'llL°°(n) 

H Vu llL~(n\s) < l|V«|U<- ( n) < 1 _i (r | +r .> ^ ■ ( 67 ) 
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Moreover, by (|64[). we have 



1 



|Vu|| L oo( B) < -llVvHioo^) 



1.1 ^C\\VS n g\\ 



< T mm {n — TT' 1 /-; — i i i ♦ ^ 

It ~~ -M 1 — M + r v e 

< C\\VS Q g\\ Lao{n) 
~ 1 — \a\ + r*^fe 

Since llViSnglU 00 ^) < C||<7||c a (dn)> we finally get ()13[). This completes the proof of 
Theorem 11.31 □ 
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